Low-frequency vibrations of soft colloidal glasses 
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We conduct experiments on two-dimensional packings of colloidal thermosensitive hydrogel parti- 
cles whose packing fraction can be tuned above the jamming transition by varying the temperature. 
By measuring displacement correlations between particles, we extract the vibrational properties of 
a corresponding "shadow" system with the same configuration and interactions, but for which the 
dynamics of the particles are undamped. The vibrational spectrum and the nature of the modes are 
very similar to those predicted for zero-temperature idealized sphere models and found in atomic 
and molecular glasses; there is a boson peak at low frequency that shifts to higher frequency as the 
system is compressed above the jamming transition. 



PACS numbers: 63.50 -x, 63.50 Lm, 82.70 Dd 

Crystalline solids are all alike in their vibrational prop- 
erties at low frequencies; every disordered solid is disor- 
dered in its own way. Disordered solids nonetheless ex- 
hibit common low-frequency vibrational properties that 
are completely unlike those of crystals, which are domi- 
nated by sound modes. Disordered atomic or molecular 
solids generically exhibit a "boson peak," where many 
more modes appear than expected for sound. The excess 
modes of the boson peak are believed to be responsible 
for the unusual behavior of the heat capacity and ther- 
mal conductivity at low-to-intermediate temperatures in 
disordered solids [2]. 

It has been proposed that a zero-temperature jamming 
transition may provide a framework for understanding 
this unexpected commonality [3]. For frictionless, ideal- 
ized spheres this jamming transition lies at the threshold 
of mechanical stability, known as the isostatic point [31 [3] ■ 
As a result of this coincidence, the vibrational behavior 
of the marginally jammed solid at densities just above 
the jamming transition is fundamentally different from 
that of ordinary elastic solids ^SrtSj . A new class of low- 
frequency vibrational modes arises because the system is 
at the threshold of mechanical stability [TUj ; these modes 
give rise to a divergent boson peak at zero frequency [5j. 
As the system is compressed beyond the jamming transi- 
tion, the boson peak shrinks in height and shifts upwards 
in frequency [5.. Generalizations of the idealized sphere 
model suggest that the boson peaks of a wide class of 
disordered solids may arise from proximity to the jam- 
ming transition [51ll2|. Moreover, the jamming scenario 
predicts that systems with larger constituents such as 
colloids should also have boson peaks. 

Colloidal glasses offer signal advantages over atomic 
or molecular disordered solids because colloids can be 
tracked by video microscopy. Vibrational behavior has 
been explored in hard-sphere colloids |13| and vibrated 



granular packings [13], but difficulties with statistics [T3] 
or micro-cracks jl4j were encountered. In contrast, we 
use deformable, thermosensitive hydrogel particles to 
tune the packing fraction in situ. Our experiments 
show unambiguously that the commonality in vibrational 
properties observed in atomic and molecular glasses ex- 
tends even to colloidal glasses, in striking confirmation 
of the jamming scenario. 

We study disordered colloidal solids composed of 
poly(A^-isopropylacrylamide) or NIPA microgel particles. 
NIPA particles swell with decreasing temperature, and 
can thus be tuned at fixed number density from a loose 
packing of small particles to a jammed packing of highly 
deformed, larger particles over the range of a few degrees 
in temperature. Particles were loaded between two glass 
cover slips, creating a monolayer in which the spheres 
were confined to move in the horizontal plane. A bi- 
nary mixture (diameters 1 fim and 1.4 /im at T = 24.7°C, 
with a large/small number ratio ^ 0.7) was employed to 
suppress crystallization. The sample was hermetically 
sealed using optical glue (Norland 63), and annealed for 
two hours at 28° C. Data were acquired using standard 
bright field video microscopy at temperatures ranging 
from 24.7°C to 27.2°C. Over this range the thermal en- 
ergy is essentially constant, so that temperature serves 
primarily to change particle diameter and packing frac- 
tion. The sample temperature was controlled by thermal 
coupling to the microscope objective (BiOptechs), and 
the sample was allowed to equilibrate for 15 minutes at 
each temperature before data acquisition. The region of 
interest was far from any boundaries or non-uniform re- 
gions within the sample. The video rate was 30 frames/s 
for a duration of 1000 s. The trajectories of the N « 3600 
particles in the experimental field of view were extracted 
using standard particle tracking techniques. Cage rear- 
rangements did not occur in any of the data sets shown 
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FIG. 1: (color online) Interaction (in terms of the effective 
spring constant kij) between large particles as extracted from 
the stiffness matrix K a.t (f> — 0.885 (squares), and (f) = 0.840 
(diamonds). The results are fit to harmonic (solid line), and 
Hertzian (dashed line) interactions, convoluted with a Gaus- 
sian size distribution to take the known polydispersity into 
account. Inset: Interaction range a and energy scale e for the 
fitted harmonic interaction as a function of packing fraction. 
The solid line in the upper inset corresponds to the particle 
diameter implied by the volume fraction. 



during the 1000 s of run time. Thus, the system remained 
in the same basin of the energy landscape and each par- 
ticle had a well-defined average position. 

As in |13l m], we define u(t) as the 2Af-component 
vector of the displacements of all particles from their av- 
erage positions, and extract the time-averaged displace- 
ment correlation matrix, or covariance matrix, as |13l 114] 



(1) 



where z, j = 1, . . . , IN run over particles and coordinate 
directions, and the average runs over time frames. 

In the harmonic approximation, the displacement cor- 
relations, C, are directly related to the stiffness matrix, 
K, defined as the matrix of second derivatives of the ef- 
fective pair interaction potential with respect to particle 
displacements. To quadratic order, the effective potential 
energy of the system is ^ = }^\\^ K\i. Within the energy 
basin the system is thermally equilibrated, so we can cal- 
culate correlation functions from the partition function, 
Z (X J Duexp (— i/Su'^iiTu), where /3 — l/k-QT. In par- 
ticular, {uiUj) — k-QT {K^^^^ 

We now introduce a shadow system of particles ther- 
mally equilibrated in the same configuration with the 
same interactions as the colloids. In contrast to the col- 
loids, which are suspended in water and whose motion 
is strongly damped, the virtual particles of the shadow 
system are undamped. The real and shadow systems 
are characterized by the same correlation and stiffness 
matrices, C and if, because these are static equilibrium 
quantities. For the shadow system, however, the stiffness 
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FIG. 2: (color online) The vibrational density of states &t <j) = 
0.859 averaged over (a) different ratios 2N/T of the number 
of degrees of freedom to the number of frames, as labeled, 
and (b) different values of the optical resolution relative to 
the average rms displacement, e/c, as labeled. 



matrix is directly related to the dynamical matrix 



Di 



(2) 



whose eigenvectors correspond to the vibrational modes. 
Thus, this analysis allows direct comparison of the 
damped colloidal solid to disordered atomic or molecular 
glasses and to idealized sphere packings. 

One advantage of this analysis is that we can obtain 
vibrational information about the disordered solid with- 
out knowing the effective interaction potential between 
particles in advance. As an added bonus, we can extract 
information about the interactions. The elements of stiff- 
ness matrix K contain direct information about the sec- 
ond derivative of the effective pair potential, denoted as 
kij. We find that kij is far above the noise only for adja- 
cent particles. This observation gives us confidence in the 
method; clearly, elastic interactions between neighboring 
particles dominate the data, as they should. Further- 
more, positional information about the large particles is 
sufficiently accurate to yield the stiffness as a function 
of separation, k{r). We fit k{r) to the second deriva- 
tive of the harmonic (a = 2) or 3D Hertzian (a = 5/2) 
interaction potentials given by 



V{nj) = 



1 - 



Ri + Ri 



<R^ + R 



(3) 



and V — otherwise. We convolute the fitting functions 
with a Gaussian distribution to account for the known 
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FIG. 3: (color online) (a) The vibrational density of states relative to the Debye prediction, Z)(aj)/a;, for (j) = 0.859. The position 
of the maximum in this plot, shown by the vertical dashed line, defines the boson peak frequency, uj* . (b) The density of states, 
D{lj) and (c) the participation ratio p{y>). The dotted vertical lines in (b) and (c) indicate the frequencies of the modes shown 
in the remaining panels, (d-f) Displacement vector plots of eigenmodes at (d) low, (e) intermediate, and (f) high frequencies. 
The size of each arrow is proportional to the displacement of the particle at that position, weighted by its mass: | u,: ] y/rrii . Blue 
circles in (d) indicate regions of high-displacements in the quasilocalized mode. 



colloidal particle polydispersity. The results are plotted 
m Fig. [1] Note that we cannot differentiate between the 
harmonic and Hertzian fits; the difference is tiny except 
for extremely large overlaps where statistics are limited. 
The fit parameters e and a = {Ri + Rj) for the har- 
monic potential are shown as insets in Fig. [T] As one 
might expect, the average size of the particles increases 
with (/), and the energy scale is roughly constant when 
4> > 0.84. The solid curve in the top inset of Fig.[l]shows 
the value of the diameter implicit in the assignment of 
(j) obtained by interpolation of dynamic light scattering 
data at different temperatures at low concentration. Here 
we have measured the elastic contact distance while dy- 
namic light scattering measures a hydrodynamic radius. 
The uncertainty in the solid curve is larger than the dif- 
ference between the curve and our data. The correlation 
method therefore provides a useful new way of measuring 
interparticle interactions at high concentrations. 

The accuracy of our results depends on the statistics 
of the time averages to calculate the matrix elements in 
Eq. [T] as well as on the ability to resolve particle dis- 
placements u(i). The optical resolution, e « 5nm, at 
best in our case, must be small compared to the av- 
erage root-mean-squared displacement, calculated to be 

1/2 

c = (Cj/ ) w 22nm at the packing fraction (p = 0.859. 
Likewise, the number of degrees of freedom, 2N k, 7200, 
should be small compared to the number of time frames. 



T = 30000, at most in our case, over which Cij is av- 
eraged. Fig. [2]ja) and (b) show the density of vibra- 
tional states, D{ll!), as e/c and 2N/T are varied. -D(w) 
is normalized so that D{uj)dLu = 2N, and e is var- 
ied by rounding measured particle displacements. Note 
that poor resolution artificially lowers the high end of 
the spectrum but that inadequate statistics raise the high 
end, so that the two effects tend to cancel. The spectrum 
appears reasonably close to convergence for T = 30000 
and e = 5nm, used in the other figures. The error in the 
width of the spectrum for these values is estimated as 
10% [m (see also [H]). 

The vibrational spectrum for crystals has the Debye 
form, D(lij) ^ lo. In Fig.jsJ^a) we plot the measured den- 
sity of states relative to the Debye prediction, D{u!)/uj, 
a.t (f) = 0.859, well above the jamming transition at 
(j)c ~ 0.84 [16] [23]. For an elastic solid such as a crystal, 
D{u!)/uj would be flat at low frequencies; the presence 
of a maximum in Fig. [3]ja) indicates the existence of a 
boson peak and the position of the maximum defines the 
boson peak frequency, uj* . Thus our results show that 
a boson peak — an important feature in the vibrational 
spectrum of atomic and molecular glasses — also appears 
in a disordered colloidal solid. 

Fig.[3]Jb) shows the vibrational density of states, D{uj). 
Fig. [3] shows the mode participation ratio, p{u}), which 
measures the degree of spatial localization of a mode n: 



4 




10 

(JO (rad/s) 



FIG. 4: (color online) (a) D{u))/u} and (b) participation ratio 
vs. frequency as a function of packing fraction (as labeled in 
(b)). Inset to (a): the boson peak frequency uj* as a function 
of packing fraction. 



P{^n) = (Ej'^i|en,jP)^/(A^Ej'^l|en,i|''), where e„_, is 
the polarization vector in mode n and is the mass, 
of particle i. Thus, p(a;„) ~ \/N for a localized mode 
and p{uJn) ~ for an extended mode. Note that the 
results are qualitatively the same even if we do not weight 
the sums with the particle masses. 

At low frequency, Fig. [sf^c) shows that the modes are 
clearly very different from plane-wave sound modes be- 
cause they have a low participation ratio; Fig. [sj^d) de- 
picts a typical low-frequency mode. It is quasi-localized 
with localized structure superimposed on a plane-wave- 
like background. At intermediate frequencies, eigen- 
modes are highly disordered and extended with p ~ 
0(1) [3 Hlj (Fig. ^e)) while at high frequencies the 
modes are localized, as expected (Fig. [sj^f)). 

Fig. Qa) shows the density of states relative to the De- 
bye prediction, D{uj)/uj, as a function of (j) and the boson 
peak frequency, w* (inset to Fig. |4|^a)). As expected, lo* 
increases with ^ It does not vanish at (pc because 
the system is at nonzero temperature. Similar thermal 
smearing of the zero-temperature jamming transition is 
observed in idealized spheres [18| . 

Fig. Qb) shows that the shape of p{uj) is similar at 
all packing fractions. As the system is decompressed to- 
wards (/)c, however, the low-w quasilocalized modes shift 
down in frequency, edging closer to instability [8]. 

Quasi-localized modes at low frequency have been ob- 
served in idealized sphere packings as well as in atomic 
or molecular glasses (see [S] and references therein). In 
the sphere packings, these modes are associated with the 



lowest energy barriers for rearrangements [S] and are the 
ones that tend to go unstable. It has been shown that 
under mechanical load, a mode shifts downwards in fre- 
quency to zero, signaling a rearrangement [12]. Such 
modes are quasi-localized, and the rearrangement oc- 
curs in the region in which the mode has high displace- 
ments (circled in Fig. [sjjd)) f5I?. Low- frequency quasi- 
localized modes are also connected to irreversible particle 
rearrangements driven by thermal fluctuations above the 
glass transition temperature [51] [52]. Thus, our experi- 
ments demonstrate that it is possible, in a real disordered 
solid for which we do not know the interparticle interac- 
tions a priori, to identify localized regions that are likely 
to be predisposed towards failure. 
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